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ABSTRACT 

We present four spectra of the Type la supernova (SN la) 2006D extending 
from —7 to +13 days with respect to B-band maximum. The spectra include 
the strongest signature of unburned material at photospheric velocities observed 
in a SN la to date. The earliest spectrum exhibits C II absorption features 
below 14,000 km s _1 , including a distinctive C II A6580 absorption feature. The 
carbon signatures dissipate as the SN approaches peak brightness. In addition to 
discussing implications of photospheric-velocity carbon for white dwarf explosion 
models, we outline some factors that may influence the frequency of its detection 
before and around peak brightness. Two effects are explored in this regard, 
including depopulation of the C II optical levels by non-LTE effects, and line-of- 
sight effects resulting from a clumpy distribution of unburned material with low 
volume-filling factor. 
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Subject headings: supernovae: general — supernovae: individual (SN 2006D) 



Introduction 



Type la supernovae (SNe la) make valuable standard candles because of their intrinsic 
brightness and the homogeneity of their light curves. Observations of high-redshift SNe la are 
responsible for the recent revelation that the rate of expansion of the Universe is accelera ting 
rtPerlmutter et alill998l : iGarnavich et alill998l : iRiess et al.lll998l : IPerlmutter et al. l ll999h . A 
push to better calibrate SNe la as distance indicators and control systematics in future high- 
precision cosmology experiments has helped motivate the search for deeper physical insight 
into SN la progenitors and their explosion mechanism. 

Certain recent multidimensional hydrodynamical SN la explosion models involve the 
thermonuclear disruption of a Chandrasekhar-mass C-0 white dwar f (WD) by deflaqra 



tion, the propagation of a subsonic flame front through the star (e.g., iGamezo et al.l 12003 



Ropke fc Hillebrandtl 120051 ) . A generic characteristic of these new models is the presence of 
unprocessed WD material below the cano nical 14,000 km s" 1 cutoff s een in W7, a tuned 



spherically-symmetric deflagration model (INomoto. Thielemann fc Yokoilll984l ). Thus, car- 
bon detection at such low velocities would generally provide observational support for the 
newer multidimensional deflagration models. Few SN la spectra actually exhibit a robust low- 
velocity carbon signature; to date the lowest measured are the C II A6580 and A 7234 features 



in th e spectra of SN 1998aq, which may extend deeper than 11,000 km s 1 ( Branch et al. 



20031 ). The lack of ubiquitous carbon signatures in the spectra of SNe la, and the observed 
deficit of kinetic energy and synthesized nickel indicate that these models are not viable, 
though increasingly sophisticated models (e.g., with greater resolu t ion an d more detailed nu- 
cleosynthesis) may prove otherwise. Most recently, iMarion et al.l (120061 ) used near-infrared 
spectroscopy to conclude that nuclear burning in at least three normal SNe la is complete be- 
low 18,000 km s~ , a conclu sion supporting models where an initial deflagration transitions 
to a supersonic detonation flKhokhlovl Il99ll ) . 



We present Nearby Supernova Factory (SNfactory. lAldering et al.ll2002l ) spectroscopy of 
the SN la 2006D, which includes an unambiguous photospheric- velocity C II signature. The 



data were obtained usin g the Supernova Integral Field Spectrograph (SNIFS. lAldering et al. 



20021 ; lLantz et al.ll2004l ) on the University of Hawaii 2.2-meter telescope on Mauna Kea. Our 
focus here is on the early photospheric-phase spectra of SN 2006D; an analysis involving all 
of our data will appear elsewhere. 
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2. Data & Analysis 



The Brazilian Sup ernova Search (IColesanti et al.ll2006l ) discovered SN 2006D in MCG 
-01-33-34 (z = 0.00853. [Davoust fc Continilbooj ) on 2006 January 11.2 UTC. Using a SNIFS 
spectrum obtained Ja nuary 14.6 UTC , we classified SN 2006D as a SN la one week prior to 
maximum brightness (lSNfactoryll2006l ) . This and three subsequent spectra appear in Figured] 
and are summarized in Table [TJ The spectra w ere reduced u s ing ou r dedicated data-reduction 
proced ure, similar to that presented in § 4 of iBacon et al.l (120011 ) and § 2 of lAldering et al. 



(120061 ) . No correction for interstellar reddening has been applied. The intermediate-mass 
element (IME) spectral features typical of premaximum SNe la are labelled on the first 
spectrum for reference. A preliminary light curve fit derived from SNIFS acquisition images 
suggests a 5-band peak date of 2006 January 21.8 UTC. The light curve is more narrow than 
is typical for a spectroscopically normal SN la — the derived time-scale "stretch" parameter 
s ~ 0.75 is significantly smaller than in the normal case (s = 1). This value is similar to 
that derived fro m the light curves of SNe 1986G and 1992bo (0.74 and 0.73 respectively, 
Guv et alJ^OOS). bu t not as extreme as in the sub- lum inous SNe 1991bg or 1999by (s ~ 0.6; 



Ruiz 



-Lapuentel 



20041 ; iGarnavich et all 12004 |jhall2006[ ). 



The dotted vertical lines in Figure [T] indicate the rest-frame wavelengths of the four 
C II lines that would be the strongest under the assumption of local thermodynamic equi- 
librium (LTE) at 10,000 K. The shaded regions to the blue of each line cover Doppler shifts 
between 14,000 and 10,000 km s -1 , representing velocities typical of the photosphere before 
and around peak brightness. A striking feature of the —7 day spectrum is the V-shaped 
absorption centered at 6320 A (rest frame), which we attribute to C II A6580 at a Doppler 
shift of 12,000 km s -1 . In the first spectrum, two much weaker notches fall within the bands 
of the shaded regions corresponding to C II AA4745 and 7234. These two features and the 
6580 A notch weaken by day +13. 

Another interesting feature in our spectra is the notch centered at 4100 A, just to the red 
of the small absorption typically attributed to Si II AA4128,4131. If the 4100 A absorption 
is due to C II A4267, then it extends from 14,000 km s _1 down to 8,000 km s -1 , slower than 
the red edge of the A6580 absorption by 2,000 km s _1 (of course, such a low velocity edge 
could be the result of line-blending). While the other C II features dissipate by day +13, 
the 4100 A notch persists and the absorption along its red edge strengthens somewhat. In 
later spectra not shown here, this feature steadily becomes washed out by iron peak lines 
in this part of the spectrum. If this feature does arise from C II A4267, then it does so 
nonthermally — in LTE at 10,000 K it would be approximately 30 times weaker than C II 
A6580; we return to this point shortly. 



In Figure[2J we plot our earliest spectrum together with premaximum spectra of SNe 1998aq, 
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1986G, and 1999by. The spectrum of SN 1998aq is an excellent match to that of SN 2006D, 
particularly redward of 4400 A, aside from the stronger C II A6580 absorption in the latter. 
Blueward of 4400 A, the agreement is not as good; the 4100 A notch in the spectrum of 
SN 2006D has no counterpart in that of SN 1998aq save for a minor wiggle (in fact, this 



wiggle has been previously observed in the spectra of SNe la, see Figure 2 of iBranch et al. 



20061 ). The spectra of SNe 1986G and 1999by exhibit larger Si II 5800 A to Si II 6150 A 
absorption ratios than does SN 2006D, indicativ e of lower temperature than in more spec- 
troscopically normal SNe la ( Nugent et al.lll995l ). The region from 4000 to 4400 A in these 
two spectra is dominated by Ti II, also characteristic of lower temperature. The differences 
between the spectra of SN 2006D and SN 1998aq in this region, and the possibility based on 
its light curve that SN 2006D is slightly sub-luminous, suggest that Ti II may play a role in 
the formation of its spectrum. 

More detailed SN spectral line identification requires an accounting for line blending, 
so we compare the —7 day spectrum of SN 2006D with spectra generated by the highly- 
parameterized SN spectrum synthesis code SYNOW (e.g.. lBranch et al.ll2003l ). SYNOW uses 
a simplified model of a SN atmosphere consisting of a sharply-defined continuum-emitting 
core surrounded by a line-forming region. Line formation is treated under the Sobolev 
approximation, with line opacity parameterized as a function of ejection velocity and with 
relative strengths for a given ion set by assuming Boltzmann-factor population of the levels. 
A synthetic spectrum is overlaid on the observed —7 day spectrum in Figure [3l 

The dashed fit reproduces most of the main absorption features of the spectrum, includ- 
ing those of IME ions Ca II, Si II, and S II. A velocity at the photosphere of 11,000 km s^ 1 
produces satisfactory results. The blend at 4700 A consists of lines of Fe II and Si II, while 
the blend from 4100 A to 4400 A consists of lines from Fe III, Mg II, and Si III, and in- 
deed some Ti II (which also fits the absorption at 3600 A as well). This Ti II detection is 
consistent with cooler, fast-declining SNe la. 

The sharpness of the A6320 feature, and the lack of a discernable accompanying emission 
feature to its red motivates us t o distribute C II Sobole v line opacity in a layer "detached" 
from the SYNOW photosphere (jJeffery fe Branch! Il990l ). An absorption feature without a 
corresponding P Cygni emission m ay also be indicative of a clump along the line of sight 
partially obscuring the photosphere (IThomas et al. A thin shell of Sobolev line opacity 

immediately above the photosphere would result in a flat absorption trough instead of a 
sharper downward spike. In the synthetic spectrum, C II opacity is detached from the 
photosphere to a velocity of 12,000 km s~ , producing a good match to the C II AA6350 
and 7234 features. The dotted synthetic spectra in the insets of Figure [3] are from a pure 
C II spectrum, overlaid on the neighborhoods of the observed C II AA4745, 6580, and 7234 
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features. 

We can account for the 4100 A feature with a separate higher velocity Ti II component 
(included in the fit with v > 17, 000 km s _1 ) but ultimately reject this hypothesis due to the 
lack of any sign of high- velocity components to explain other spectral features. In particular, 
Ca II opacity should trace the same high- velocity processed material, and we detect no such 
signature. 

In the absence of any obvious alternative to Ti II for the 4100 A notch, we suggest 
that the feature could come from C II A4267, stimulated by some non-LTE process. To 
investigate this hypothesis, we examined the non-LTE departure coefficients (6, = rii/nf TE , 
where rii is the population number of the atomic l evel, i) of the lower levels of the optical C II 



lines in a PHOENIX calculation of the model W7 (IBaron et al.ll2006l ). In the C-O-rich zone, 
departure coefficients in optical lines are found to be significantly smaller than unity (10~ 7 
to 10~ 3 ). Furthermore, the ratio of the A4267 lower-level bi to that of A6580 is observed to 
range between 1 and 60. This situation arises due to more efficient recombination through 
the ultraviolet resonance line C II A687, so that recombination into the optical levels is 
suppressed and they take on level populations significantly smaller than LTE. The size of 
this ratio could account for the persistence of C II A4267 while the other C II lines dissipate. 
Further detailed models are required to discover whether this effect is indeed generic and 
holds in other models besides W7. This effect may also resolve the long-standing difficulty 
of making a reliable identification of C II in the spectra of SNe la. 



Discussion 



That some SN la spectra possess C II signatures, particularly before maximum light, 
i s now quite c lear. Past tentati ye identifications of C II A6580 in such cases as SN 1990N 



(IMazzalil l2001[ ) and SN 1999ac flGaravini et al. 



2005 J, and previous ly unidentified notches 



seen in SN 1994D and SN 1996X pointed out by lBranch et al.l (120031 ) now seem more mean- 
ingful. These detections and their circumstances prompt consideration of several issues 
relevant to the study of SN la explosion models. 

The primary issue is the amount of carbon detected. To place a lower limit on the 
carbon mass (particularly in light of the small frj's observed in the PHOENIX model), we 
compute LTE electron level populations and corresponding C II A6580 Sobolev optical depth 
on a temperature-density grid, assuming a C-O-rich composition. The absolute minimum 
density giving rise to optical depth r ~ 1 at 10 days after outburst is 5.9 x 10~ 15 g cm -3 . 
Treating this as an average density results in a lower limit of 0.007 M of carbon (0.014 M 
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of unburned material) over the velocity range 10,000 to 14,000 km s 1 where the absorption 
feature is detected. This value is smaller than the 0.085 M & of unbur ned material presen t 



at the same velocity interval in the multidimensional model bl50_3d of iRopke et al.l (120061 ). 



If the behavior of the C II b^s seen in the C-O-rich zone of the W7 PHOENIX model is 
generic, then the amount of unburned material detected could be significantly higher than 
this estimate. The fej's need not be so miniscule as observed in the PHOENIX model — here 
only a factor of 6 is needed to produce general agreement with model bl50-3d. Extremely 
small departure coefficients drive the mass up but the mass is obviously limited by the mass- 
energy budget of the progenitor, fixed by the Chandrasekhar mass. One clear way to make 
the detection consistent with this constraint is for the unburned material to be distributed 
in clumps if the fej's are indeed very small. 

The presence of unburned material at the observed velocities may appear more consistent 
with published multidimensional deflagration models than with delayed-d etonation models, 



but the situation is not quite so simple. For example, iHoflich et al.l (120021 ) demonstrate that 
sub-luminous delayed-detonation models can fit the light curve and infrared spectra of the 
sub-luminous SN 1999by, and account for observed C I features observed with a small amount 
of carbon extending below 14,000 km s _1 . Interestingly, the best-fit model for SN 1999by is 
from a grid of delayed-detonation models with variable transition density, and corresponds to 
the lowest luminosity /lowest transition density in the grid. The trend in those models seems 
to indicate that less luminous models, with lower transition densities, have a correspondingly 
decreased carbon cutoff velocity. Thus, if the more luminous SN 2006D were to fall on that 
particular grid of models, it would need to correspond to a higher transition density and 
have a higher carbon cutoff velocity than SN 1999by, at odds with our observations. More 



speculatively, iGamezo. Khokhlov fc Oranl (120041 ) outline a scenario in which a small amount 



of unburned material in pockets could be left by a multidimensional delayed detonation 
- as the detonation wave processes material left between plumes generated during the 
deflagration phase, abrupt twists and turns could cut it off from fresh WD fuel. Elucidating 
the details (the mass and spatial distribution of such pockets) will rely on progress in future 
multidimensional delayed-detonation models. 

It is interesting to consider the significance of the sporadic and diverse nature of the 
detections of the C II A6580 feature itself. The premaximum rest-frame 5800-6600 A spectra 
of SN 2006D and five other SNe la plotted in Figure H] exhibit a variety of strengths and 
velocity ranges in C II A6580. The theoretical multidimensional models themselves may 
suggest an explanation for this phenomenon beyond any non-LTE effect — the distribution of 
unburned material at these velocities could be quite asymmetric, confined to clumps of small 
filling factor compared to the photosphere. Thus, differences from SN to SN would naturally 
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be explained as a function of the distribution of the clumps with respect to the observer's 



l ine of sight. This argument is consistent with the complementary finding by iThomas et al. 



( 120021 ). where the relatively consistent shape and depth of Si II A6355 absorption at maximum 
light indicates that the distribution of IME's in SNe la cannot be highly asymmetric. Of 
course, differences in the details of temperature structure and ionization from one SN to 
another may be in play, though we would expect temperature effects to run along a continuum 
correlated to, say, radioactive nickel mass. If low-velocity carbon detections are sporadic and 
less correl ated with luminosi ty (note that the light curve of SN 1998aq is characterized as 
"typical," iBranch et al.ll2003l ). one would expect the geometrical argument to be at least as 
plausible. An underlying mechanism for either explanation (or both) will depend on further 
observations of the earliest SN la spectra that probe the outer layers of these events, and on 
improved detailed modeling of these spectra. 



4. Conclusions 

We have presented the strongest evidence to date of unburned ejecta at low velocity 
in the early optical spectra of a SN la. The case of SN 2006D reiterates the importance 
of early-time spectroscopy (and spectropolarimetry) for SN la research. Complimentary 
observations in the nebular phase (one year after explosion) may present further opportunities 
to probe the distribution and amount of unburned material at low velocity in SN 2006D; 
while their predictions are not strictly applicable at the velocit i es con sidered in this work, the 



multidimensional spectrum synthesis models of iKozma et al.l (120051 ) suggest that unburned 
material at velocities below 10,000 km s _1 roughly consistent with our mass estimates may 
be detectable at late times as forbidden [O I] and [C I] emission. 
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Table I. Journal of Spectroscopy of SN 2006D 



Julian Day 


UTC Date 


Exposure 


Airmass 


Median 


S/N a 


Conditions 




(2006) 


(s) 




blue 


red 




2453750.10 


Jan. 14.60 


900 


1.20 


76 


85 


non-photometric 


2453752.12 


Jan. 16.62 


1000 


1.16 


72 


69 


cloudy 


2453756.07 


Jan. 20.57 


1800 


1.22 


205 


136 


near Moon 


2453770.11 


Feb. 03.61 


1000 


1.16 


72 


93 


photometric 



a Blue channel: 3500-5050 A with 2.4 A bins. Red channel: 5000-9000 A with 3.0 A bins. 
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Fig. 1. — SNIFS spectroscopy of SN 2006D. Phases are expressed relative to a 5-band peak 
brightness date of 2006 January 21.8 UTC. Dotted vertical lines mark rest wavelengths of 
C II lines AA4267, 4745, 6580, and 7234. Dark bands indicate blueshifts between 10,000 and 
14,000 km s _1 with respect to these lines, typical of the velocity at the photosphere at these 
phases. The inset is a zoom of the region around the 6320 A notch in the —7 day spectrum. 
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SN 1998aq -9d 
SN 2006D -7d 
SN 1986G -5d 



SN 1999by -4d 



4000 5000 6000 7000 

Rest Frame Wavelength (A) 

Fig. 2. — Comparison of the —7 day spectrum of SN 2006D to premaximum spectra of 
SNe 1998aq, 1986G, and 1999by. 
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Fig. 3. — Synthetic spectra overlaid on the —7 day spectrum of SN 2006D. Inset plots show 
the coincidence between the synthetic C II spectrum and the data. 
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Fig. 4. — Comparison of C II A6580 feature across several premaximum SN la spectra. The 
spectra are resampled into 10 A bins. 



